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SUMMARY 

Parametric studies were performed to assess the effects of various parameters on the free-vibration behavior 
(natural frequencies) of [±0]2 angle-ply, fiber composite, thin shell structures in a hygrothermal environment. 
Knowledge of the natural frequencies of structures is important in considering their response to various kinds of 
excitation, especially when structures and force systems are complex and when excitations are not periodic. The 
three-dimensional, finite element structural analysis computer code CSTEM was used in the Cray YMP 
computer environment. The fiber composite shell was assumed to be cylindrical and made from T300 graphite 
fibers embedded in an intermediate-modulus, high-strength matrix. The following parameters were investigated: 
the length and the laminate thickness of the shell, the fiber orientation, the fiber volume fraction, the tempera- 
ture profile through the thickness of the laminate, and laminates with different ply thicknesses. The results 
indicate that the fiber orientation and the length of the laminated shell had significant effects on the natural 
frequencies. The fiber volume fraction, the laminate thickness, and the temperature profile through the shell 
thickness had weak effects on the natural frequencies. Finally, the laminates with different ply thicknesses had 
an insignificant influence on the behavior of the vibrated laminated shell. Also, a single through-the-thickness, 
eight-node, three-dimensional composite finite element analysis appears to be sufficient for investigating the 
free-vibration behavior of thin, composite, angle-ply shell structures. 


INTRODUCTION 

High-speed flight vehicles require new materials able to withstand hostile environments (moisture and 
temperature) at speeds that are a multiple of the speed of sound and improved computational approaches to 
designing them. The new materials, known as elevated-temperature composites, have high strength, are light in 
weight, and can be tailored for the required performance. These materials are already finding applications and 
acceptance in aircraft frame and engine structures. Cost-effective usage of elevated-temperature composite 
materials requires new ideas and innovative concepts in material selection and structure fabrication in order to 
design them for durability by properly accounting for the many failure mechanisms inherent in composites. The 
influence of hostile environments, nonlinear material, and structural behavior and the coupling between 
responses induced by various loads require sophisticated analysis methods. 

A stand-alone, multidisciplinary computer code, CSTEM (Coupled Structural, Thermal, and Electromagnetic 
Analysis and Tailoring) (refs. 1 to 7), has been developed by integrating three-dimensional, finite element, 
structural and stress analysis methods with several single-discipline codes including those for integrated 
composite mechanics (refs. 8 and 9). The advantage of performing finite element analysis using three- 
dimensional brick elements is that the temperature gradient as well as the stress gradient can be varied through 
the thickness of the laminated structure, in contrast to two-dimensional plate elements. 

This report presents parametric studies that used the CSTEM computer code to assess the effects of various 
parameters on an aircraft engine composite shell structure subjected to free vibration at elevated temperatures. 
Knowledge of the natural frequencies of structures is important in considering their response to various kinds of 
excitation, especially when structures and force systems are complex and when excitations are not periodic 



(refs. 10 and 11). The structure studied is a thin, cylindrical composite shell. The composite is to be made from 
[±0] 2 angle-ply laminates containing plies with uniform thickness. Parametric studies were performed to 
e xamine the effects on the natural frequencies of parameters, such as the length and laminate thickness of the 
shell, the temperature profile through the thickness of the laminate, the fiber volume fraction, the fiber 
orientation, and laminates with different ply thicknesses. 


BRIEF DESCRIPTION OF SIMULATION PROCEDURE 

The general-purpose computer code CSTEM (ref. 4) computationally simulates the coupled multidisciplinary 
structural, heat transfer, vibration, acoustic, and electromagnetic behavior of elevated-temperature, layered, 
multimaterial composite structures. All the disciplines are coupled for nonlinear geometrical, material, loading, 
and environmental effects. CSTEM is based on a modular structure (fig. 1) and is accessed by the user through 
its executive module. The structural response can be predicted at all the composite scales including constituent 
(fiber and matrix) and ply level. 

The composite mechanics is simulated with the ICAN (Integrated Composite Analyzer) module (refs. 8 
and 9), which performs, among others, through-the-thickness point stress analysis. The micromechanics 
equations embedded in ICAN include the effects of temperature and moisture. ICAN simulates the behavior of 
polymer composites from the constituent material level to the laminate level, as depicted in figures 2 and 3. 
ICAN includes a resident material property data bank for commercially available fiber and matrix constituent 
materials at room temperature. The user needs only to specify a code name for the desired constituent materials 
(rather than having to manually input all the properties) in the CSTEM input 

The procedure for computationally simulating composite structures is shown in figure 2 and consists of four 
parts: the constituents, the synthesis, the finite element structural analysis, and the decomposition. For a detailed 
description of these four parts see references 8 and 9. The computational procedure steps are as follows: 

(1) Constituents . The material properties of the matrix in any operational temperature and moisture 
conditions are updated by using the constitutive equation depicted in figure 3. The fiber properties were 
assumed to be unchanged by the hygrothermal conditions. 

(2) Synthesis . The thermal and mechanical properties of the plies are synthesized from the mechanical 
properties of the fibers and matrix by using composite micromechanics theory. The extensional stiffness matrix 
A, the coupling stiffness matrix B, and the bending stiffness matrix D are computed by using laminate theory 
and the thermal and mechanical effective laminate engineering properties (ref. 12). 

(3) Finite element analysis . The finite element analysis is performed at the composite structural level, and 
am ong others, the resulting forces and moments are computed at each node of the finite element mesh. Because 
ICAN performs through-the-thickness point stress analysis, steps (1) and (2) and the upcoming step (4) are 
performed at each individual node of the finite element mesh. 

(4) Decomposition . The stresses and strains in each ply are computed by using laminate theory. The stresses 
and strains within the matrix and the fibers for each ply are computed by using composite micromechanics. 

The CSTEM code is written in FORTRAN programming language and is installed on the NASA Lewis 
Cray YMP computer system. 
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SUMMARY OF STRUCTURAL EQUATIONS 


The equations of equilibrium governing the linear dynamic response of the structure written in matrix form 
are the following (refs. 10 and 11): 


[M]U + [D]U + [K]U = {R{t)} 


( 1 ) 


where 

[M] structural mass matrix 

[D] structural damping matrix 

[K] structural stiffness matrix 

{/?(/)} external applied load vector, as a function of time 
{ U } displacement vector 

{(/} velocity vector 

{ 0} acceleration vector 


The equilibrium equations (1), in the case of the free- vibration problem, are simplified by setting the load vector 
equal to zero, {/?(/)} = 0. If the damping effect is not taken into account, the damping matrix is equal to zero. 


[D] = 0. 


The equilibrium equations are simplified as follows: 

[M]U + [K ]U = 0 


( 2 ) 


The structural matrices [M] and [K] are computed by discretizing the continuous structure into a number of 
finite elements. If [H] is the shape function of a finite element of the structure, then 

(1) The mass matrix of the element is given by 

[m] = Jp [H} T {H)dV 


where p is the mass density of the element. 

(2) The stiffness matrix of the element is given by 

[k] = f[B] T [C][B]dV 

V 


where [C] is the material matrix of the element, which contains the anisotropic effective properties of the 
laminate and is related to the stress-strain relationship by 

{O} = [C] {e} 


The matrix [B] is related to the strain displacement relationship by 
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{£} = [B] {«} 


where [B] = [L] [N], [L] is a linear operator, and [N] is the matrix that consists of the shape functions [H] 

(ref. 13). 

The structural mass matrix [M] is assembled from the element mass matrices [m ]j for; = 1, n el , where n d is 
the number of finite elements in the structure. Similarly, the structural stiffness matrix [K] is assembled from 
the element stiffness matrices [k ]j,j = 1, n tV 

After the computation of the [M] and [K] matrices the next step is to solve equation (2). The solution of 
equation (2) can be postulated to be of the following form: 

U = {<{>} sin to (r - f 0 ) (3) 


Substituting equation (3) into equation (2) gives 

[K] {<>} = 00 2 [M] {*} 


(4) 


For the computation of the vibration frequencies (eigenvalues) and mode shapes (eigenvectors), CSTEM 
incorporates the following two methods: the determinant search and the subspace iteration (refs. 1 1 and 14). The 
search method is intended more for a small number of equations; the subspace iteration is more efficient for a 
large number of equations. 


GEOMETRY, LOAD HISTORY, AND FINITE ELEMENT MODEL OF LAMINATED SHELLS 

The thin composite shell had a cylindrical geometiy with the following dimensions: ratio of inner radius to 
lamina thickness of the shell Rlt = 33.3; R = 25.4 cm; ratio of longitudinal length to radius UR = 2, 4, and 6. 
The boundary conditions were one end fixed and the other free (fig. 4). The laminate consisted of continuous 
fibers maH* of the graphite material T300 embedded in an intermediate-modulus, high-strength (IMHS) matrix. 
The material properties of the fiber and matrix were taken from the data bank available in CSTEM (tables 1 
and 2). The fiber volume fraction (FVR) was 55 percent and the moisture content was 2 percent. A balanced 
and nonsymmetric laminate was assumed with fiber orientations of [±6]2 and plies of equal thickness. 

The simulation of the processing of the laminated cylinder was taken into account. The glass transition 
(curing) temperature T of the matrix was initially 215.55 °C (table 2) and was gradually reduced until the 
operating temperature was equal to 7^/2 (fig. 4). At the end of the processing, residual stresses were induced in 
the constituent materials of the composite structure. 

For the free-vibration analysis of the shell structure the effect of damping was neglected. The effect of 
interlaminar shear stresses was taken into account and is described in detail in the CSTEM users manual (ref. 4). 
The computational simulation was performed by using the three-dimensional, finite element analysis code 
CSTEM. One element (eight nodes) through the thickness of the composite shell structure was used. The finite 
element mesh consisted of 760 nodes and 360 elements (fig. 5). 
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RESULTS AND DISCUSSION 


In this section the results obtained for the different [±0] 2 angle-ply composite shells are presented and 
discussed. The parameters investigated include 

(1) Effects of fiber orientation 

(2) Effects of shell length 

(3) Effects of laminate thickness 

(4) Effects of temperature 

(5) Effects of fiber volume fraction 

(6) Effects of different ply thicknesses 


Effects of Shell Length 

The influence of shell length in conjunction with ply angle 0 was examined to study the free-vibration 
behavior of the composite shell. Results obtained for natural frequencies and mode shapes are summarized here. 

Natural frequencies . — Three thin composite shells with length ratios UR of 2, 4, and 6 were examined. The 
computed values of the first natural frequency are plotted versus 0 and the three UR ratios in figure 6. The com- 
puted values are normalized with respect to the maximum frequency of 554.9 Hz. Correspondingly, the 
computed values of the second natural frequency are normalized with respect to the maximum frequency of 
609.6 Hz and are plotted in figure 7. Finally, the computed values of the third natural frequency are normalized 
with respect to the maximum frequency of 871.8 Hz and are shown in figure 8. The following conclusions were 
reached: 

(1) For all UR the natural frequencies increased for 0° < 0 < 22.5° and decreased for 22.5° < 0 < 90°. The 
maximum frequencies (first, second, and third) occurred at 0 = 22.5°. 

(2) The shorter the shell, the higher the frequencies, as expected because shorter shells are stiffer than 
longer shells. 

(3) The percentage differences of the natural frequencies between the short (UR = 2) and long (UR = 6) 
shells for the different ply angles are significant (table 3) and depend on the ply angle. 

Mode Shapes . — The behaviors of the mode shapes of the composite shell with UR = 6 and RJt = 33.3 for 
the first three frequencies are presented graphically as follows: 

(1) First (lowest) frequency: The mode shape (fig. 9) is primarily bending, indicating that the shell behaves 
lik e a cantilever beam. 

(2) Second frequency: The mode shape (fig. 10) is a combination of torsional and breathing (or radial) 
modes. 

(3) Third frequency: The mode shape (fig. 11) is a combination of torsional, breathing (or radial), and axial 
(or longitudinal) modes. 

These mode shapes are summarized in table 4 for 0° < 0 < 90°. 
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Effects of Laminate Thickness 


The influence of laminate thickness of the thin composite shell versus the ply angle 0 for laminate thickness 
ratios RJt of 20, 33.3, and 100 is presented in figure 12. The computed values are normalized with respect to the 
maximum frequency of 264.4 Hz. The following conclusions were reached from figure 12: 

(1) For all RJt the first natural frequencies increased for 0° < 0 < 15° and decreased for 15° < 0 < 90°. The 
maximum frequency occurred at 0 = 15°. 

(2) The thicker the composite shell, the higher the values of the natural frequencies. Thus, thick composite 
shells (RJt = 20) had the highest frequencies, and thin composite shells (RJt = 100) had the lowest frequencies. 
The frequencies maximized at 0 = 15°. 

(3) The percentage difference of the natural frequencies between the thick laminated shell (RJt = 20) and the 
thin laminated shell (RJt = 100) was relatively small (table 5). 

Therefore, the laminate thickness ratio RJt had only a small effect on the first natural frequencies of the 
laminated shell. 


Effects of Temperature 

The effects of temperature and moisture profiles through the laminate thickness on the free vibration in 
conjunction with ply angle 0 were examined. Three temperature profiles with uniform moisture (2 percent) were 
considered. The first profile was uniform with temperature equal to room temperature (21 °C), the second 
profile was also uniform with temperature equal to TJ2 (107.7 °C), and the third profile varified linearly, with 
temperature equal to T ^2 at the inner surface and room temperature at the outer surface of the composite shell. 
The length ratio UR was 4 and the laminate thickness ratio RJt was 33.3. The computed values of the first 
natur al frequency, normalized with respect to the maximum value of 276.8 Hz, are plotted versus ply angle 0 
for three temperature profiles in figure 13. The following observations were made: 

(1) For all temperature profiles the natural frequencies increased for 0° < 0 < 20° and decreased for 
20° < 0 < 90°. The frequency peaked at 0 = 20°. 

(2) The room-temperature composite shell had the highest frequencies, followed by the composite shell with 
the linear profile and finally the composite shell with the uniform temperature equal to TJ2. 

(3) The percentage difference in natural frequency between the uniform room-temperature profiles and TJ2 
was low (table 5) and was almost independent of 0. 

Therefore, the temperature profile had little effect on the natural frequencies of the composite shell throughout 
the angle-ply range. 


Effects of Fiber Volume Fraction 

The effects of fiber volume fraction in conjunction with ply angle 0 on the free-vibration behavior of the 
composite shell were evaluated. Three fiber volume fractions were examined: 55, 60, and 65 percent. The length 
ratio UR was 4 and the laminate thickness ratio RJt was 33.3. The computed values of the first natural 
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frequency, normalized with respect to the maximum value of 278.6 Hz, are plotted versus ply angle 0 for the 
various fiber volume fractions in figure 14. The following conclusions were reached: 

(1) For all fiber volume fractions the natural frequencies increased for 0° < 0 < 20° and decreased for 
20° < 0 < 90°. The maximum natural frequency occurred at 0 = 20°, which was the same for all the cases 
investigated. 

(2) The composite shell with 65-percent fiber volume fraction had higher frequencies than the shell with the 
60-percent fiber volume fraction, which had higher frequencies than the shell with 55-percent fiber volume 
fraction. 

(3) The percentage difference in the natural frequencies between the two laminates with fiber volume 
fractions of 65 and 55 was very low (table 5). 

Therefore, the fiber volume fraction had little effect on the natural frequencies of the composite shell. 


Effects of Ply Thickness 

The effects of laminates with various ply thicknesses in conjunction with ply angle 0 on the natural 
frequencies of the composite shell were examined. The length ratio LJR was 4, the laminate thickness ratio RJt 
was 33.3, and the fiber volume fraction was 55 percent. Four laminated shells with different ply thicknesses and 
the same laminate thickness t were examined: (1) a [±0] angle-ply laminate with equal ply thicknesses (f/2, r/2), 
(2) a [±0] 2 angle-ply laminat e with nonuniform ply thicknesses (3f/8, f/8, t! 8, 3f/8), (3) a [±0] 2 angle-ply 
laminate with uniform ply thicknesses (r/4, r/4, r/4, r/4), and (4) a [±0] 4 angle-ply laminate with equal ply 
thicknesses (r/8, r/8,...,r/8). The computed values for the first natural frequencies of these examined laminates are 
plotted versus 0 in figure 15. The following conclusions were reached: 


(1) The natural frequencies increased for 0° < 0 < 20°, decreased for 20° ^ 6 < 90°, and peaked at 0 — 20°. 


(2) The laminas with different ply thicknesses had an insignificant effect on the natural frequencies of the 
thin composite shells. 


FINITE ELEMENT ANALYSIS USING EIGHT-NODE, THREE-DIMENSIONAL 
ELEMENTS VERSUS FOUR-NODE PLATE ELEMENTS 

The behaviors of the shell structure subjected to free vibration using eight-node, three-dimensional elements 
versus four-node plate elements were compared. The geometry of the laminated structure consisted of a cylinder 
with length equal to 40 in., an inner radius equal to 10 in., and uniform laminate thickness equal to 0.3 in. The 
laminate consisted of graphite T300 fibers and a high-strength IMHS matrix. The fiber volume fraction was 
55 percent, the ply angle was [0/90] 2 , and the moisture was 2 percent. The laminated structure had an initial 
temperature of 215.55 °C that was gradually reduced to the operational temperature of 108 °C. The natural 
frequencies were computed by using the MHOST computer code (ref. 15) in conjunction with plate element 75, 
normalized, and plotted in figure 16. The frequencies obtained with the eight-node, three-dimensional elements 
were higher than the frequencies obtained with the four-node plate elements. The reason is that the eight-node 
element is stiffer than the plate element and therefore the natural frequencies are higher. 
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In older to improve the accuracy of the computational simulation of laminated, thin shell structures 
subjected to free vibration, it is recommended that a 20-node brick element be used or that an eight-node brick 
element be used and the mesh refined. 


SUMMARY OF RESULTS 

The free vibration of thin composite shell structures in a hygrothermal environment was computationally 
simulated by using the CSTEM three-dimensional, finite element analysis computer code. The simulation 
included the effects of parameters such as ply angle, shell length and laminate thickness, temperature, and fiber 
volume fraction on the free-vibration behavior of the composite shells. The important results are summarized as 
follows: 

(1) The ply angle and the length of the laminated shell significantly influenced the free-vibration (frequen- 
cies and mode shapes) behavior of the composite shells. 

(2) The fiber volume fraction, the laminate thickness, and the temperature profile through the laminate 
thickness had little effect on the free-vibration behavior of the composite shells. 

(3) The laminates with different ply thicknesses had an insignificant effect on the free-vibration behavior of 
the composite shells. 

(4) A single through-the-thickness, eight-node, three-dimensional-composite finite element analysis appears 
to be sufficient for investigating the free-vibration behavior of thin, composite, angle-ply cylindrical shells. 
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TABLE 1.— T300 GRAPHITE FIBER PROPERTIES 
AT ROOM TEMPERATURE 


Modulus in longitudinal direction, GPa 220.7 

Modulus in transverse direction, GPa 13.79 

In-plane Poisson’s ratio 0.20 

Out-of-plane Poisson’s ratio 0.25 

In-plane shear modulus, GPa 8.96 

Out-of-plane shear modulus, GPa 4.827 

Thermal expansion coefficients in longitudinal 

direction, 10“** mm/(mm-°C) 0.99 

Thermal expansion coefficients in transverse 

direction, 10 -6 mm/(mm-°C) 10.00 

Thermal conductivity in longitudinal 

direction, W/(m-K) 83.69 

Thermal conductivity in transverse 

direction, W/(m-K) 8.369 

Fiber tensile strength, GPa 2.413 

Fiber compressive strength, GPa 2.069 


TABLE 2. — IMHS MATRIX PROPERTIES 
AT ROOM TEMPERATURE 


Modulus in longitudinal direction, GPa 3.448 

Modulus in transverse direction, GPa 3.448 

In-plane Poisson’s ratio 0.35 

Out-of-plane Poisson’s ratio 0.35 

Thermal expansion coefficients in longitudinal 

direction, 10 -6 mm/(mm-°C) 64.8 

Thermal expansion coefficients in transverse 

direction, 10 -6 mm/(mm-°C) 64.8 

Thermal conductivity in longitudinal 

direction, W/(m-K) 2.16 

Thermal conductivity in transverse 

direction, W/(m-K) 2.16 

Matrix tensile strength, MPa 103.4 

Matrix compressive strength, MPa 241.3 

Matrix shear strength, MPa 89.6 

Void fraction 0.225 

Glass transition temperature, °C 215.55 


TABLE 3 — PERCENT DIFFERENCE* IN NATURAL FREQUENCIES BETWEEN 
SHORT (UR = 2) AND LONG (UR = 6) CYLINDERS 


Frequency 

[0] 4 

[±15] 2 

[±30] 2 

[±45] 2 

[±60] 2 

[±75] 2 

[±90] 2 

Difference in natural freqencies, percent 

First 

57 

60.7 

63.5 

70.7 

73.3 


75.4 

Second 

60 

59.6 

67 

65.9 

63.5 


40 

Third 

58 

45.4 

50 

55.4 

54.9 

l 

57.3 


“Percent difference = = 6) * 100%. 
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TABLE 4.— MODE SHAPES OF LAMINATED SHELL 
{UR = 6, /e/r = 33.3, AND FVR = 55%) 

AS A FUNCTION OF PLY ANGLE 


Ply angle, 

e, 

deg 

First 

frequency 

Second 

frequency 

Third 

frequency 

f±e ] 2 

Beam mode 

Torsional and 
radial modes 

Torsional, axial, and 
radial modes 


TABLE 5.— PERCENT DIFFERENCE IN FIRST NATURAL FREQUENCIES 


Laminates 

(0] 4 

[±15] 2 

[±30] 2 

[±45] 2 

[±60] 2 

[±75] 2 

[±90] 2 



Difference in natural freqencies, percent 


Thick!#// = 20) and 
thin {RJt = 100) a 

2.35 

2.5 

4.5 

9.62 

1.8 

0.8 

1.2 

Room-temperature 
{T = 21 °C) and 
high-temperature 
(7= T/Z) b 

3.8 

2.5 

3.1 

3.6 

3.3 

3 

3.5 

FVR of 65% and 
FVR of 55% c 

6.06 

6.05 

6.1 

6.3 

6.2 

5.4 

5.8 


a Percent difference = 


Frequency!#/? = 20) - Frequency!#/? = 100) innC w 
frequency!#// = 20) 


Percent difference = 


Frequency! T = 21 °C) - Frequency (T ~ TJ2) 

„ — — x 100% 

Frequency (T ~ 21 °Cj 


Percent difference = 


Frequency (FVR = 65%) - Frequency(FVR = 55%) v ]()()<%, 
Frequency (FVR = 65%) 



Figure 1 . — CSTEM modular structure. 
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Figure 2. — Integrated Composite Analyzer (ICAN). 
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Figure 3. — Regions of constituent materials and 
nonlinear material characterization model. 
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Figure 4. — Geometry, materials, and processing of 
laminated thin shells, (a) Geometry: UR = 4; R/t = 33.3; 

R = 25.4 cm. (b) Materials: angle-ply laminate [0/-8] 2 ; 
fiber volume ratio = 55%; moisture = 2%; T300 fibers; 
IMHS matrix, (c) Processing. Thermal load T g = 215.5 °C. 
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Normalized second frequency 




Figure 5. — Three-dimensional finite element mesh and transverse 
cross-sectional area of laminated thin shells. 



Figure 6. — Influence of shell length ratio on lowest frequency 
of laminated shell. R/t = 33.3; FVR = 55%; moisture = 2%; 
maximum normalized frequency = 554.9 Hz. 



Figure 7. — Influence of shell length on second frequency of lamina- 
ted shell. R/t = 33.3; FVR = 55%; moisture = 2%; maximum 
normalized frequency = 609.6 Hz. 



Figure 8. — Influence of shell length on third frequency of laminated 
shell. R/t = 33.3; FVR = 55%; moisture = 2%; maximum normal- 
ized frequency = 871 .8 Hz. 


13 



Figure 9. — First mode shape of [±0]2 angle-ply laminated shell. 


Figure 1 0. — Second mode shape of [±0]2 angle-ply laminated shell. 



Figure 1 1 .—Third mode shape of [±0]2 angle-ply laminated shell. 


Figure 12. — Influence of laminate thickness ratio R/t on lowest 
frequency of laminated shell. UR = 4; FVR = 55%; moisture = 2%; 
maximum normalized frequency = 264.4 Hz. 
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0 10 20 30 40 50 60 70 80 90 

Ply angle, 6, deg 

Figure 13. — Influence of temperature profile on lowest frequency of 
laminated shell. UR = 4; R/t = 33.3; FVR = 55%; moisture = 2%; 
maximum normalized frequency = 276.8 Hz. 



Ply angle, 6, deg 


Figure 14. — Influence of fiber volume fraction on lowest frequency 
of laminated shell. UR = 4; R/t = 33.3; moisture = 2%; maximum 
normalized frequency = 278.6 Hz. 



0 10 20 30 40 50 60 70 80 90 

Ply angle, 0, deg 

Figure 1 5. — Influence of laminate thickness on lowest frequency of 
laminated shell. UR = 4; R/t = 33.3; FVR = 55%; moisture = 2%; 
maximum normalized frequency = 262.3 Hz. 



Mode 


Figure 16. — Comparison of finite element analyses for eight-node, 
three-dimensional elements and four-node plate elements. 
Maximum normalized frequency = 454.1 6 Hz. 
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